Quantum Transport in Ladder- Type Networks: Role of nonlinearity, topology and spin 
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We investigate quantum transport of electrons, phase solitons, etc. through mesoscopic networks 
of zero-dimensional quantum dots. Straight and circular ladders are chosen as networks with each 
coupled with three semi-infinite leads (with one incoming and the other two outgoing). Two trans- 
mission probabilities (TPs) as a function of the incident energy e show a transition from anti-phase 
aperiodic to degenerate periodic spectra at the critical energy Ec which is determined by a bifurca- 
tion point of the bulk energy dispersions. TPs of the circular ladder depend only on the parity of 
the winding number. Introduction of a single missing bond (MB) or missing step doubles the period 
of the periodic spectra at e > Ec . Shift of the MB by lattice constant results in a striking switching 
effect at £ < £c- In the presence of the electric-field induced spin-orbit interaction (SOI), an obvious 
spin filtering occurs against the spin-unpolarized injection. Against the spin-polarized injection, on 
the other hand, the spin transport shows spin-fiip (magnetization reversal) oscillations with respect 
to SOI. We also show a role of soliton in the context of its transport through the ladder networks. 

PACS numbers; 03.75.-b, 05.45.-a,05.60.Gg. 



I. INTRODUCTION 

Recently there has been a growing interest in quan- 
tum transport in discrete physical systems characterized 
by networks with nontrivial topologies P, Q . Those net- 
works mimic networks of nonlinear waveguides and and 
optical fibers Q , Bose-Einstein condensates in optical 
lattices 0], superconducting ladders of Josephson junc- 
tions [B], double helix of DNA, etc. In these networks, 
their topology and the presence of a few embedded de- 
fects are expected to play a vital role in controlling the 
macroscopic quantum transport such as a switching of 
the network current. Here, a main interest lies in the net- 
works connecting everywhere-discrete lattice points 0] 
in contrast to another topical works on quantum graphs 
which are composed of connected continuous linear seg- 
ments of finite length ■ 

On the other hand, with introduction of the nonlin- 
earity to the time-dependent Schrodinger equation, the 
network provides a nice playground where solitons prop- 
agate in a complicated way until escaping through the 
attached semi-infinite leads. There already exists an ac- 
cumulation of studies of the soliton propagation through 
the discrete chain, and its collision with small defect clus- 
ters 0. However, little work has been done on the soliton 
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transport through the big networks with and without de- 
fects. 

In this paper we investigate quantum transport of elec- 
trons or phase solitons through mesoscopic networks of 
zero-dinicnsional quantum dots. Typically, straight and 
circular ladders are chosen as model networks with each 
being coupled with three semi-infinite leads (with one 
incoming and the other two outgoing). In Section II , 
based on the discrete cubic nonlinear Schrodinger equa- 
tion, we examine a fate of the soliton coming from the 
incoming lead and propagating through the above net- 
works in a complicated way until escaping through the 
three semi-infinite leads. The two transmission probabil- 
ities (TPs) based on a soliton picture are evaluated and 
compared with the result of Landauer formula based on 
the (stationary and discrete) linear Schrodinger equation. 
The following Sections are based on the standard (linear) 
quantum mechanics. In Section III, TPs are explored as 
a function of the incident energy, and the characteristic 
features of the transmission spectra are found. In Section 
IV we shall elucidate a radical change of the transmission 
spectra by introducing a single defect bond into the net- 
work. The role of topology in the transport through the 
circular ladder is also studied in this Section. Finally in 
Section V the electric-field induced spin-orbit interaction 
(i.e., Rashba interaction) is introduced to the network. 
Then we investigate the result of spin transport through 
the networks and indicate its role in magnetization oscil- 
lations and spin filtering. Summary and discussion are 
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devoted to Section VI. 

II. MODEL NETWORKS AND DISCRETE 
NONLINEAR SCHRODINGER EQUATION 

As a challenge to analyze general big networks, we 
choose two type of networks, straight and circular lad- 
ders (see Figs. [T]and[2]), which mimic Josephson junc- 
tion or double helix of DNA. Each system consists of 
an array of zero-dimensional quantum dots (i.e., lattice 
sites), where central part represents a network and exter- 
nal three lines stand for the attached semi-infinite leads. 
All lattice points are numbered in the way given in Figs. 
[Hand [21 In Fig. [l] for example, the incoming lead (left) 
is connected with the ladder at the site m and a pair 
of outgoing leads (right) are connected with it at the 
sites m -f 2n and m + 2n + 1 . Suppressing three external 
leads, the ladder includes 2n + 2 lattice sites and n — 1 
steps (perpendicular to the ladder). The wave function 
comes through the incoming lead (^in), collides with the 
network, and is partly reflected through the incoming 
lead (*&re/) and partly transmitted through two outgo- 
ing leads {^outi,^out2)- Dynamics of a wave function 
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FIG. 1: Straight ladder with 3 leads. 




FIG. 2: Circular ladder with 3 leads, 
in these open networks is described by discrete nonlinear 



Schrodinger equation (DNLSE), 



(1) 



where A represents the strength of cubic nonlincarity. Aij 
is adjacency matrix giving the topology of the network 
and is defined, in a suitable energy unit (say, K) by 



if j and I are linked 
otherwise 



(2) 



In the case of quantum dots with a common discrete level 
(CDL) for each, (t) is the wave function of the j-th dot. 
The distances between linked lattice sites are fixed to a 
common value, say, d with d being of order of 10 ~ 100 
nm. K stands for the tunneling matrix element between 
connected adjacent dots. CDL is chosen around Fermi 
energy and prescribed to zero energy. Time t is in units 
of h/2K and A = U/2K with U the very weak Hartree 
term due to the electron-electron interaction. Firstly we 
investigate the injection of a wave packet (WP) through 
the incoming lead, where DNLSE governs: 



(3) 



Consider, at i = 0, Gaussian WP centered at ^o, with 
initial momentum kg and width 70 . In its discrete version 
the time-dependent WP can be written as 



$j (t) = y/N exp 



1' 



+ ^k{j - + ^7;ij - 0' 



(4) 

where ^(i) and j(t), which are scaled by d, are time- 
dependent center of mass and width of WP, respectively. 
k{t) and 5{t), which are scaled by d"^ and d~^, respec- 
tively, are the corresponding canonical-conjugate vari- 
ables. 

In the limit •jd ^ d, WP dynamics can be obtained 
from effective Lagrangian 



i = fee - 7'- - 




cos(/c)e 



(5) 



from which we have the equations of motion for ^, fc, 7 and 
6. In order to have a stable WP (soliton) on incoming 
leads it should be 7 = 5 0, from which it follows 



A. 



cosfc 
70 



(6) 



with ^ < k{= ko) < TT and (5 = 0. Under this condi- 
tions we present the numerical results of soliton dynamics 
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FIG. 3: Soliton dynamics in straight ladder with 3 leads. 
Time evolution of the spatial distribution of the positive wave- 
function probability: 1^2^ (3,3') (4,4'). k = fyr. 
Basal Lengths and wave number are scaled by d and , 
respectively. Ladder steps are not depicted for simplicity. 
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FIG. 4: Soliton dynamics in circular ladder with 3 leads. 
Time evolution of the spatial distribution of the positive wave- 
function probability: 1-^2-^ (3,3') (4,4'). k = fyr. The 
same notion on lengths, wave number and ladder steps holds 
as in Fig|3] 



colliding with a network in Figs. [3]and[H Soliton propa- 
gates through the incoming lead (marked as '1'), collides 
with network (marked as '2'), propagates through net- 
work (marked as '3' and '3") and is partially reflected 
through the incoming lead (marked as '4') and partially 
transmitted through two outgoing leads (marked as '4" ) . 
Transmission and reflection probabilities (TP and RP) 
at long enough time after collision with the network can 



be calculated as 
R 



E 

outgoing lead 1 

E 

outgoing lead 2 

E 



(7) 
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The result as a function of the incident wave number 
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FIG. 5: Comparison Ti,T2 and R between Eq. Q with use 
of nonlinear dynamics of a soliton and Eq. in Landauer 
formula for the time-independent linear Schrodinger equation. 
Number of steps in ladder is n = 10. Solid line and '-f ' for 
Ti, dashed line and 'x' for T2, and dotted line and for R. 

k (scaled by d"^) is shown in a set of symbols in Fig. [5] 
in the case of the straight ladder with number of steps 
n — 10 and length of each external lead m = 250. Here 
initial width of wave packet 70 = 50 and initial center of 
mass Co = 100. We find the unitarity T1+T2 + R=l 
is always satisfied, namely no fraction of WP remains in 
the central network at long-enough time. 

Also, we compare this result with the result based 
on Landauer formula applied to the time- 

independent linear Schrodinger equation for the ladder 
network with N{= 2n) lattice sites, which is connected 
with the semi-infinite incoming lead at '0' site and two 
semi-infinite outgoing leads at 'N + V and W -I- 2' sites. 
In the latter approach, the outgoing wavefunction 
($0, • ■ ■ , <^N+i,'^N+2 f is determined by [T^ 



(8) 



against the incoming wave function ^'j„ = 
{-K,[F-\+) - F-i(-)]$o(+),0,...,0)^ with Ks 
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and F the tunneling and transfer matrices, re- 

spectively, in the leads. G is the Green function defined 

by 



G 



E- H 



(9) 



III. TRANSMISSION SPECTRA OF STRAIGHT 
LADDER 

One cannot recognize any universal feature in Fig. O 
in the case of a ladder with ri = 10 steps. However, when 
n 3> 10, there appear universal characteristic features 
independent of n. In Fig. [6] transmission and reflection 



In Eq.®, H is the Hamiltonian which includes the in 



teraction of the network with external leads 
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where 
Vo,Vn 



H is 

^1, Vn+2 





Vn+2 
(10) 

Hamiltonian. 



the unperturbed 
and Ko^i,Kn-i,n+i,Kn.n+2 are 
respectively the self-energies which renormalise the 
effect of semi-infinite leads and the tunneling matrices 
between the ladder network and leads. Noting that all 
tunneling matrices are unity by scaling in the present 
calculation, we reach the transmission Tj with j = 1,2 
and reflection probabilities R, 

T, ^ \< N + j\G\0 > k:[F-\+) 



a = 1,2), 




R 



\< o\G\o > k;[f-\+) 

-F-H-)]-lf 



(11) 



In Fig. Owe compare the results of Eq. ([7]) with those of 
Eq. (fTTj) in case of the ladder with N — 20. Surprisingly 
two approaches give the identical results. The reason is 
that the width of the WP employed here is much longer 
than the linear dimension of the network and that the 
nonlinearity plays little role. Precisely speaking, so far 
as the soliton is large enough and fast enough to guar- 
antee that the time of collision between the soliton and 
ladders is much shorter than the soliton dispersion time, 
one may resort to a linear approximation to compute the 
transmission coefficients [l3| . In the following, therefore, 
we shall derive Ti, T2 and R with use of Eq. pT]) applied 
to the linear Schrodinger equation for the latter. 
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FIG. 6: Transmission and reflection probabilities against in- 
cident energy e. Number of steps in the ladder is n = 50, 100 
and 200 from top to bottom panels. Solid, dashed and dotted 
lines correspond to Ti , T2 and R, respectively. Ti and T2 are 
degenerate for e > 0.5. 



probabilities against energy (e) of the incoming electron 
are plotted in case of the straight ladder with n = 50, 100 
and 200 steps. The unitarity Ti + T2 + i? = 1 is always 
satisfied. We find the existence of a critical energy Ec = 
0.5 and the remarkable difference of TPs between the 
lower (0 < e < £c) and higher {sc < s < 1) energy 
regions. In the lower energy side, Ti and T2 have the 
anti-phase structure (i.e., Ti takes peaks whenever T2 has 
dips and vice versa), and the oscillation period decreases 



5 



as £ — > Ec- In the high energy side, on the other hand, 
two TPs are degenerate and highly periodic. All these 
characteristics hold irrespective of the value of n, so long 
as the network is big enough (n ^ 10). In fact, we 
obtained the same spectrum in case of n = 1000 as in 
Fig. [6l while the oscillation period is further shortened 
in the latter. 

The mechanism underlying the above characteristics 
is explained by using the perturbation theory. Let's first 
investigate the nature of the unperturbed long network 
without three leads, which can be regarded as a periodic 
ladder in Fig. [7l For a pair of upper and lower sites 2m 



2m+1 



-cos (X) -0.5^— 
-cos (X) to .5 




FIG. 7: Unperturbed periodic straight ladder. 



and 2m + 1, the wave functions satisfy 

£^•2™ = -2(*2(m+l) + *2(m-l) + *2m+l), 
£$2m+l = -^(*2(m+l) + l + *2(m-l) + l + «'2m).(12) 

Let us introduce new basis functions Um and Vm with use 
of the transformation: 



Um = -J^{^2m + <&2m+l) 
Vn 



7^^ 

75(^2™ - $2m+l) 



(13) 



Um and w„i stand for the even- and odd-parity states 
in each step, respectively. Using this new basis, the 
eigenvalue problem is decoupled, namely, reduced to the 



even- and odd-parity parts. Assuming Ur, 



and 



ues 



for an infinitely long ladder, we find eigenval- 



Eu = - COs(fc) - - 

E^ = -cos(fc) -I- i. 



(14) 



The even-parity branch e„ and odd-parity one Ey con- 
stitute a pair of energy bands (see Fig. [8|). It should be 
noted: while for < e < Ec, both energy branches and 
Ey appear, only the Ey branch can survive for e > Ec- 




FIG. 8: Two branches of energy dispersion for unperturbed 
ladder. Vertical axis stands for energy e. £„: even-parity 
branch; £„: odd-parity branch. 



Under the presence of the perturbation, namely, in the 
case of the ladder attached with three leads in Fig. [2 
Uni,Vm ^ e**^™ are not the eigenstates any more: the 
mixing (superposition) of states occur within the odd- 
parity manifold only for e > Ec and between the odd- 
and even-parity manifolds for < e < Ec- In case of e > 
Ec, therefore, the wave function retains the same feature 
as the unperturbed state: the coefficients of the wave 
function $2m and $2m+i have the identical magnitude. 
This fact holds at the ladder edge with m = 2n and 
TO = 2n -|- 1 as well. Consequently, we see the degeneracy 
of oscillations for Ti and T2 in Fig. [S] On the other hand, 
in case of < e < Ec, we see the superposition of Um and 

aUm + Pvm = ■^{a + f3)(^2rn + ^{a- fi)^2m.+l- (15) 

As a result, wherever the coefficient of $2m has a big 
magnitude, that oi<^2m+i has a small one, and vice versa. 
This is true even at the ladder edge, explaining the anti- 
phase oscillation for Ti and T2 in Fig. [HI 

Thus, the transmission spectra of the straight ladder 
attached with three leads show a mixing between differ- 
ent parity states and anti-phase structure in the output 
in the lower energy regime (0 < e < Ec), while, in the 
higher energy regime {Ec < £ < 1), no mixing and the 
degenerate periodic structure in the output. 



IV. ROLE OF DEFECT BONDS AND 
TOPOLOGY 

One of the most essential question of quantum net- 
works is whether or not only a single defect bond in- 
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troduced into big networks will plays a crucial role in 
quantum transport. Now we proceed to investigate the 
influence of a missing bond embedded in the midst of 
the ladder network with iV = 100 steps on the quantum 
transport. The left and right panels in Fig. [5] correspond 
to breaking a bond and step, which are parallel and per- 
pendicular to the ladder, respectively. The correspond- 
ing transmission spectra are given in Figs. [10] and [TlJ 

Missing bond Missing step 
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FIG. 9: Missing bonds {A,B) and missing steps {C,D). Each 
figure shows only 6 quantum dots in the midst of the long 
regular ladder. A{B) corresponds the case that a single bond 
with X , which is parallel to the ladder, is missing. Missing 
bond in (B) is displaced upwards from one in [A) by lattice 
constant; C{D) corresponds the case that a single step with 
X, which is perpendicular to the ladder, is missing. Missing 
step in (D) is displaced to right from one in (C) by lattice 
constant. 



Consider the case with a missing bond (MB) in the mid- 
ladder. For e > Ec, the regular oscillation of Ti and T2 
retains the degeneracy and in-phase structure, but has a 
period twice as large as the one without MB. For £ < Ec, 
Tl shows a radical change from the complete transmis- 
sion (Tl = 1) to the complete reflection (Ti — 0) and vice 
versa when MB moves by lattice constant, which can be 
taken as a switching effect (see Fig. [T0|). The issue of a 
missing step (MS) in the midst of the ladder is as follows: 
for £ > £c, besides the period-doubling phenomenon, the 
regular oscillation shows a phase shift by half a period 
when MS moves by lattice constant (see Fig. [TT|) . We 
should note: so long as a reference MB or MS is embedd 
in the midst of big networks, the above discoveries (i.e., 
period doubling and phase shift for e > Ec, and switching 
effect for e < £c) remains unchanged, irrespective of the 
absolute location of such a defect bond in Fig. O Thus, 




FIG. 10: Transmission probabilities against incident energy e 
in case of a single missing bond (MB). Numbers of steps (n) 
and of lattice points in the ladder are 100 and 200, respec- 
tively. Upper panel includes Ti (solid line) and T2 (dashed line) 
in the case that MB lies between lattice points 100 and 102 
(: case 'A' in Fig. [9}. Lower panel includes only Ti, and solid 
and dashed lines correspond to cases 'A' and 'B' in Fig. [9] 
respectively. Spectra are degenerate for e > 0.5. 



an introduction of a single MB or MS into a big network 
results in a radical change in the transmission spectra. 

In order to see the role of another topology of networks 
we consider the annular circular ladder and investigate 
the twist effect (see Fig. [T^ on quantum transport. 

In the case of no twist, the spectra show the same re- 
markable transition when s crosses £c = 0.5 as in the case 
of the straight ladder. We find: In the lower energy side, 
Tl and T2 have the anti-phase structure, and the oscilla- 
tion period decreases as £ ^ Sc- In the high energy side, 
on the other hand, two TPs are degenerate and highly pe- 
riodic. In the presence of a single twist (i.e., analogue of 
Mobius strip) the spectra again shows a remarkable tran- 
sition at £c = 0.5, but the detailed feature differs from 
the result for the no twist case. See the great reduction 
of Tl and T2 in the lower energy region in the single twist 
case. On the other hand, in the double twists case the 
result is identical to that of no twist case. The spectra is 
determined by the parity of the winding number (WN). 
The winding of the circular ladder is identical to the ap- 
plication of Aharonov-Bohm flux with WN multiplied by 
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FIG. 11: The same as Fig. [10] but in the case of a single 
missing step (MS). Upper panel includes Ti (solid line) and 
T2 (dashed line) in the case that MS lies between lattice points 
100 and 101 (: case 'C in Fig. (Qjl. Lower panel includes only 
Ti, and solid and dashed lines correspond to cases 'C and 'D' 
in Fig. [9] respectively. 




FIG. 12: Twisted annular circle. ' x ' means the disconnection, 
h + 1 and h + 2 (likewise, h and h + 3) axe connected. 

a half of the flux quantum ^ = ff • Thus the topology 
of networks plays a vital role in quantum transport. 



V. SPIN-ORBIT INTERACTION AND SPIN 
TRANSPORT 

Recent progress in semiconductor spintronics revealed 
a way of controlling the magnetization of devices not by 
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(a) No twist 
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(b) Single twist 
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FIG. 13: Transmission probabilities against incident energy e 
(solid for Ti and dashed for T2). Three cases of twisted circles: 
a) no twist; b) a single twist; c) double twists. Spectra are 
degenerate for e > Sc, though no bold line appears. 



a magnetic but by an electric field. The_idea is to use 



14 



15 



Rashba spin-orbit interaction (SOI) 
strength is tuned by the external gate voltage. In this 
Section, by introducing SOI into the network, we inves- 
tigate spin transport (spin-dependent transport) as well 
as charge transport. According to the pioneering work 
of Datta and Das [ll|, [3i 13 , we first consider the spin 
transport against the spin-polarized injection. The net- 
work Hamiltonian generalized so as to include Rashba 
SOI is given by 



-^Aj-i$,-fa(CTXp)^$j=e$j (16) 



with <I>j = {(f'j.h 't'j-l)'^ the two component wave func- 
tion, a = — ^^f^-ij^ Ez the strength of Rashba SOI in the 
case of an vertically applied electric field and a stands 



for Pauli matrices. In Eq. (|T6|) . energy is scaled by the 
tunneling matrix element K. For convenience in our 
numerical calculation, we introduced dual ladders to as- 
sign each of them to up- and down-spin states, respec- 
tively (see Fig. [T4l) . The spin transport is quantified as 



spin 



^1,2(1) ^ ^1,2(1) and the charge transport as 



In Fig. [15] the spin transport against incident energy 
is plotted for different values of the strength of Rashba 
spin-orbit interaction a. We consider the spin-polarized 
[Sz — +\) injection. In the absence of spin-orbital inter- 
action the spin transport (STP) 7^^*", T^^*" as a function 
of £ show the same spectra as in the case of charge trans- 
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FIG. 14: Spinor ladders. For computational purpose, dual 
ladders are introduced with each corresponding to spin-up 
and spin-down states. 

port Ti, T2 (see Fig. |6]), because we have no contribution 
from Ti^2(i)- Against the variation of SOI, the spin trans- 
port shows spin-flip (magnetization reversal) oscillations 
(see Fig. [15]), while keeping the anti-phase structure of 
rpspm rpspm range e < £c(= 0.5). Against the 

variation of SOI, by contrast, the charge transport (CTP) 
keeps the spectral feature without SOI (see Fig. [S]). 

Finally we shall investigate the most interesting sub- 
ject, namely the spin transport in network systems with 
SOI against the injection of spin-unpolarized electron. 
Figure [1^] shows T^p'^'^ and T^^™ as a function of e for 
non-zero values of a. Astonishingly we find T^p^" = 
_rpspm value of £ in the case of a 7^ 0. This dis- 

covery indicates that a straight ladder with three leads 
plays a role of the spin filtering, i.e., the unpolarized elec- 
tron is decomposed into mostly spin-up and mostly spin- 
down components through its transport in the ladder. In 
the context of nanoscience, this is the most essential issue 
among many other discoveries in the present work. 



VI. SUMMARY AND DISCUSSIONS 

Choosing straight and circular ladders as big net- 
work models and attaching them with one incoming and 
two outgoing semi-infinite leads, we examined quantum 
transport of an electron or phase soliton. In the begin- 
ning, by adding a small cubic nonlinearity (e.g., Hartree 
term) to the discrete time-dependent linear Schrodinger 
equation, we showed how the incoming soliton bifurcates 
at the entrance of the ladder-type network and is ulti- 
mately evacuated from the network through three leads. 

We chose a soliton large enough and fast enough to 
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FIG. 15: Spin transport T^^^"^ (solid) and T^^'" (dashed) 
for different values of spin-orbital interaction in case of spin- 
polarized injection. The panels from top to bottom corre- 
spond to a = 0.0,0.12 and 0.18, respectively. 



guarantee the time of collision between the soliton and 
ladders to be much shorter than the soliton dispersion 
time. On the basis of this soliton picture, two trans- 
mission probabilities (2^1,2) and a reflection probabil- 
ity (i?) were evaluated, which proved to accord with 
the corresponding probabilities obtained from the linear 
methodology, i.e., Landauer formula applied to the time- 
independent linear Schrodinger equation. The main part 
of the paper was then devoted to the results of the latter 
(linear) methodology. Firstly we investigated Ti, T2 as a 
fimction of energy e of the incident electron. Both proba- 
bilities show a transition from anti-phase aperiodic to de- 
generate periodic spectra at the critical energy Sc = 0.5, 
whose value is determined by a bifurcation point of the 
bulk energy dispersions. TPs of the circular ladder de- 
pend only on the parity of the winding number (WN), 
because WN plays a role of Aharonov-Bohm flux with 
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FIG. 16: The same spin transport T^'"'"' (solid) and T^™ 
(dashed) as in Fig. 1151 but in the case of spin-unpolarized 
injection. The panels from top to bottom correspond to a = 
0.1,0.12 and 0.46, respectively. 



its magnitude being a half of flux quantum multiplied by 
WN. 



Introduction of a single defect bond into big networks 
radically changes the macroscopic transport spectra. A 
missing bond (MB) parallel to the ladder in the network 
doubles period of the periodic spectra for e > Sc- For 
e < Ec, shift of a single MB by lattice constant results 
in the switching between two outgoing leads. A missing 
step leads to a phase shift besides the period doubling 
for e > Ec- 

Finally, by introducing the electric-field-induced 
Rashba spin-orbit interaction (SOI), we explored spin 
transport (y^'"?'"^ rpspm-^ against the spin-polarized in- 
jection. At zero SOI, 7^''™ and J"*^™ as a function of e 
show the same spectra as in the case of charge transport. 
Against a variation of SOI, however, this structure shows 
a coherent spin-flip (magnetization reversal) oscillations. 
On the other hand, the injection of the spin-unpolarized 
electron leads to the spin filtering, namely, the unpolar- 
ized electron is decomposed spatially into mostly spin-up 
and mostly spin-down components through its transport 
in the ladder. Therefore the present network can be used 
as a spin-filtering device. This is the most striking issue 
of this paper. The present results would also be applica- 
ble to propagation of a wide-enough soliton in Josephson 
junction networks and of a wave packet in Bose-Einstein 
condensates in optical-lattice networks, although the lin- 
ear and static approximation will break down and the 
transport would be highly nonlinear and more generic. 
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